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ABSTRACT 

This dissertation deals with the energy method for 
the torsional analysis of prismatic members which esdiibit 
a linear or non— linear stress -strain behaviour. Various 
stress— strain relations in nondimens ional fortn are used 
as idealisation of material properties so that these can 
be used to cover wide range of concretes and other stru- 
ctural materials. 

A stress function satisfying the boundairy condit- 
ions and having unknown coefficients is assuiued^Che unkn- 
own coefficients in the stress function are determined 
by minimising the complementary potential. Results obtained 
numerically/ are presented for rectangular and circular 
cross-section for different parameters m the generalised 
Raroberg-Osgood relation, thereby covering materials like 
concrete, steel, aluminium, magnesium and reinforced 
plastics. Comparison of results with the existing experi- 
mental and analytical results is made. Finally a methodo- 
logy for ultimate strength theory using hyperbolic sine 
relation for stress strain expression is presented. 
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: area of cross-section 

: radius of circular cross-sections: 

: dimensions of rectangularcroross-sectionj, 
width = 2a., depth = 2b, 

A^^: Unknown Coefficients 

; modulus of elasticity inofeension 

: modulus of elasticity in tension for plastic 
stirain. 

; Strain 

; Assumed form of function for strain 
: cylinder compressive strength of concrete 
: tensile strength of concrete 
j Shear modules 

; Shear modules for plastic portion 
: 2a 
: b/a 

T 

m 

0 mm I I ' Mm ■ 

Sa^bf; 
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j : polar moment of interia of test specimen 

i/j/k,r/p,q: 1,2,3 

L : Length of prismatic member 

r I one of the parameters in Ramberg-Osgood function 

m : r+1/2 

s : arc parameter 

R i In the small eighbourhood of the solution 
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: Torque 
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’■y-' : Ultimate Torque 

u,v,w : displacement components 

Uj, : total complementary energy in the member 

: total complementary energy qjer unit -volume 
X/y,2 : coordinates 

x^y : x/a, y/b 
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: Total cdinplementary worlc done by the external 
torque . 

: caaplementary potential energy 
: stress function 
: twist per unit length 
: Ultimate twist per unit length 
; warping function 
: shear strain exponents 
: shear strain resultant 

: shear strain resultant for elastic and plastic 
portion respectively 

# 

: Ultimate shear strain of material 
; shear stress resultant shear stress resultant 
: shear stress resultant 

; non— dimensionalising shear stress and shear 
strain 

: parameter of Ramberg-Osgoad function 
; maximum normal tensile strain and tenisilo stress 
: Poisson's ratio 
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1. INTRODUCTION 

Methods exist for the analysis of torsion of 
prismatic members of isotropic, homogeneous materials which 
are either linearly elastic or rigid perfectly plastic. 
Structural materials in general, are neither elastic nor 
perfectly plastic throughout the range of behaviour either 
in pure compression or pure tension. This simple fact under- 
lines the necessity and importance of analytical work for 
non-linear elastic materials in torsion. 

In the past such a work was not considered feasible 
because of the difficulties the material non-linearity 
introduces. It makesthe resulting simultaneous equations 
or the governing differential equations non-linear in 
nature and at that time these non-linear equations were 
rather difficult to solve. So, naturally fca" design work 
etc., resort was taken to the then existing results by 
assuming the material to be either elastic or perfectly 
plastic. But now, with the advent of electronic computers 
along with the vast development of numerical techniques, 
hurdles created by non-linearity has altogether been 
eliminated. Though, inspite of all this, little effort has 
been made in this field of inelastic analysis for torsion 
but definitely research is gaining momentum in this field 
The present work is a step in this direction. 
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1,1 Present Method and Objectives 

An attempt is made herein to develop a non— linear 
elastic analysis for prismatic members in torsion for short 
term loading using energy approach^ for general idealisation 
of inelastic material behaviour. The objectives of an inela- 
stic analytical method are two-fold: 

1. To get an analytical undexstanding of the behaviour 
of prismatic members subject to torsion> from begining to 
failure . 

2, The knowledge of inelastic behaviour helps in the 
understanding of shear stress-strain behaviour from 
uniaxial tension, compression and torsion tests. 

In this analytical investigation various stress 
strains curves are used, such as: parabolic and polynomial 
expressions, Ramberg-Osgood function, Sinh curve. Any of 
these curve can be used to represent the shear stress - 
strain behaviour of the given material. A complementary 
energy approach is used wherein the stress function ^ is 
assumed to satisfy the boundary conditions. The unknown 
parameters in the stress function are deteemined by 
minimising the total comolementary energy? potential for 
an assumed value of the twist 0 . Then the twisting 
moment and the maximum resulting shear stress r for this 
value of 0 are determined. Thus theoritical plots of 
T Vs Q and ^ cuives are obtained. These are 

compared with the available analytical and experimental 


results . 
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1»2 Review of Previous Work 

Coulomb^ was first to give an exact solution of 
the torsional problem for a circular bar df homo^eheous, 
isotropic, elastic material* The tnairi assumption behind 


such a solution was that plane section" normal to the long- 
itudinal axis of twist always remain plane* Later when 

(2) 

Navier tried to extend the above theory to noncircular 
sections he ended with contradiction with the boundary con- 
ditions. However, such an attempt made it iclearcithat plane 
sections normal to the longitudinal axis of twist in case 
of nonculular cross-sections does not remain plane after 
the application of tongue . Saint Venant^^^ guided by Navier ‘s 
conclusion assumed both rotation and warping of section and 
put forwarded the so called semi -inverse method for the 
solution of torsion problem for prismaticu members* After- 
wards mathematical solutions for the torsion problem for 
different cross-sections have been obtained. Seth^'^^ gave 
the solution of T-bean with infinite flanges by using 
complex variable approach whereby a T-shaped section was 
mapped on a unit half circle by using the Swatz Chritoftel 
transformation. 

(5) 

lyeangar provided the rigrous solution for T,L 
and I sections. By dividing the flanged sections into two 
rectangles i.e, web rectangle and flange rectangle, he 


*Number in brackets designate references at the end. 
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assumed two stress functions for these two rectangles . 

These stress function were determined in such a way that 
besides satisfying the boundary conditions oh the respe- 
ctive rectangle, these were ,req;uired to satisfy the conti- 
nuity conditions at the junction. 

For some of the sections, the mathematical analysis 
became so complicated that it gave rise to a new method of 
solving torsion problem. This was to find mathematical 
analogy between the torsion problem and the behaviour of 
some other physical problem which is easy to visualize and 
solve. One such analogy was given by Prandtl^^^ and is 
known as Membrame Analogy. This is based on-hhe mathemati- 
cal similarity between the torsion problem and the behavi- 
our of a uniformly stretched elastic membrane subjected 
to uniform lateral pressure. The membrane has the same 
shape as the cross-section and its edges are restrained. 

The deflection of such a membrane give rise to a differ- 
ential equation similar to the governing differential 
equation in torsion problem. And hence the deflection 
function in the above problem corresponds to the stress 
function in the original problem. Such an analogy is useful 
in the e3<perimental determination of stress function and 
stress distribution. 

Based on the membrane analogy, soap films as 
membranes have also been in use for the direct measure- 
ment of stresses. Two holes are made in a plate, one of 
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the same cross-section as under investigation and the 
other a circular hole of predetermined diameter. Films 
on both these holes are applied with the sanle lateral 
pressure. By measuring the slopes of two films, comparison 
of the stresses of the given member with the circular bar 
for the same unit twist can be made. The ratio of the 
volume under the two films give the corresponding torque 
ratio . 

Besides this, there are number ofh hydrodynamical 
(7) 

analogies. Kelvin sUQ'gested the §iroilarity between stress 
function and stream function of the irrotational motion of 
an ideal fluid in a vessel of the same cross-section as 
the twisted bar, 

( 8 ) 

BousSinesq showed the similarity between deter- 
mination of stress function and Velocities in a laminar 
motion of viscous fluid in a tube of the same cross-section. 
All these hydrodynamical analogies may not be important for 
the actual solution of torsion problem, but these defini- 
tely help one to draw cettain important conclustins which 
otherwise are difficult to draw^ One such conclusioh is 
that at the projecting corners of a cross-section of a 
twisted bar the shearing stresses become zero whereas 
at the reentrant comer these stresses become theoritically 
infinitely large. 

Buchanan for the first time studied the torsion 
problem of prismatic members of isotropic, homogeneous 
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material with non-linear stresses tram relation. This 
non— linear relation between shearing stress and shearing 
strain was determined experimentally f romitoffsion tests 
on circular beams. Such a method of obtaining shear stress 
strain curve gave good results as confirmed by the compar- 
ison of computed and experimental results. 

Two series of tests were conducted. One was to 
study the failure mechanism and the cracking pattern of 
rectangular and T-shaped cross— section for various 
combinations of longitudinal and transverse reinforcement 
From these tests it was concluded that failure of rein- 
forced T beams depends upon the type and location of 
reinforcement . 

In the second series of test, three circular, 
one rectangular and one T-shaped Cross-section without 
reinforcement were loaded to failure in pure torsion to 
obtain torque-rotation dat^k The theoritical analysis 
employed is nothing but the modification of Saint-Venants 
theory to take into account the material nonlinearily • 
4|hear modulas G, which is constant for a "linear case 
now becomes a function of the stress level in otherwords 
a function of the coordinates (x,y) in the- cross-section. 
Finite difference technique was employed to solve the 
governing differential equation. Each cross-section was 
divided into small grids and a difference equation was 
written for each grid point. By solving the simultaneous 
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9 

equations, {JT was obtained for a particular value of 9. 

The 9 IS then increased by a small amount and the corres- 
ponding value of shear stress and torque ^were determined. 
Comparison with experimental results shows good correlation. 

In all these test, model beams of plaster were used 
in place of structural doncrete prototypes. This was done 
with the view to reduce the curing time but care was taken 
to take a proper mix of plaster so as to represent the 
stress strain behaviour of concrete prototype. 

From the above study, Buchanan concluded that for 
structural concrete which does not have plastic plafeau 
on its stress strain curve, strain hardendd distribution 
IS more suited than plastic stress distribution near 
failure. According to him, there is very little difference 
in the canputed torque capacity for either ultimate stress 
distribution, which possibly explain why plastic stress 
distribution have given satisfactory results in the 
past . 

Smith and Sidebottom^^*^^ made use of principle of 
minimum complementary energy potential for solving the 
torsion problem for a material with non-linear stress strain 
law. Shear stress-strain relation was represented by hyper- 
bolic Sinh function in terms of normal (tensile) stress 
strains. This was done by making use of the vena Mises 
Octahedral shear stress theory to convert uniaxial state 
of stress by multiaxial state of stress. The stress function 
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52 f satisfying the boundary conditions was obtained by 
minimising the complementary potential w.r.t. each of 
the unknown parameters in the stress function. Knowing 
shear stress and torqu4 for different value of twist 
were obtained. However such a solution is approximate 
unless the number of functions satisfying the boundary 
conditions in the stress function tends t:o infinity. 

Taking few such functions sufficiently accurate result 
were obtained for rectangular and trianguar cross-sections. 
For circular case, exact results have been_ given. 

Shah, Setlur & Chatterjee^ have applied the 
finite element method. The continuous two dimensional 
problem was discretized by dividing the cross-section 
into a finite number of triangular elements . Complementary 
potential energy of the member expressed ±n terms of the 
nodal stress function values and angle of twist, was 
minimised and thereby obtaining a set of non-linear 
algebraic equations . These were solved nuerically to give 
5? for one value of unit twist. Thus plotso of torque 
Vs twist for circular, rectangular and L cross-section 
for different material properties were obtained. Compa- 
rison of results with existing results shew good agreement. 

Design of R.C.C. structures for torsion is still 
considered a difficult problem in structures of ordinary 
proportion. With the increase in use of monolithic struc- 
tures, it has now become an important problem So lot of 
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research is Being done to obtain tobsional strength of 
concrete membters/ which is essentially a problem concerning 
non-linear material. The main probleh\ with concrete is that 
it is weak ih tension than in shear With the result that it 
fails in tension much before a failure due to shear takes 
place. The appfoximate relation betweeti compressive, ten- 
sile and sheaf strength of concrete is in the ratio of 
10:1:2 . Besides this other difficulties associated with 

it &T9 , determii>ation of maximum tensile .-stress and its 
variation and ptoblem of arbitrary cross *^sectioti. Perhaps 
these are some of the points which justify the eeprate 
treatment that hdS been given to it besides other struct- 
ural materials. 

Plain cotldfete is rarely used where torque is 
likely to be applied. But tests, discussed later, show 
that addition of reinforcement does not change the behav- 
iour of members under torsion, atleast up to the first 
crack. Knowledge of the behaviour of the section upto first 
crack is also important from practical considerations like, 
serviceability, deflection and rotation. Moreover the 

increase in torsional strength due to reinforcement is 

( 13) 

hardly of the order of 15% . Tests by Hsu showed 

that rectangular sections of plain concrete fails as soon 
as the first crack appears. Torque and twist curve show 
a linear relation upto approximately 8C% of the torsional 
capacity before begining to show slight non-linear behavi- 
the inelastic redistribution of stress near 


our due to 
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failure. 

In flanged sections of plain concrete in pure 
torsion the failure is not so sudden. Zia*" observed 
that first crack appears in the web and the beam continues 
to carry more torque until the crack propagates through 
flanges and the manber fails. 

Hsu on the basis of ej<periments concluded that 
rectangular concrete beams in pure torsion actually fail 
by bending . This was based on the criteria that failure 
is reached when the tensile stresses induced by a 45^ 
bending component on the under face of the rectangle 
reaches a reduced modules of rupture . Making use of 
equilibrium at failure^ he gave a mathematical relation 
for ultimate torque, valid for rectangular cross-section. 
For other flanged rections he suggested the algebraic 
sum of torque for each of the divided rectangle. The torque 
so given is on conservattve side and has been supported by 

experimental results . 

( 14 ) 

Zia found that basically the torsion theories 
for plain concrete can be expressed in a basic form as 

KAD 

where X = maximum Corsional stress 
T = torque 

A = Cross-sectional area 
D =: diameter of the inscribed circle 
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The summary of various theories in the tabular 

\ 

form is given in Table i values of K depend on the 

ratio of h/b and varies between 0,2 to 0.5. (—) the 

4A 

shape factor, where P represents perimeter of the section, 
And'b^ and h^^ are width and depth, respectively, of the 
component rectangles . 

Research on R.C.C. members subjected to torsion 

by number of authours< 13- 14-20) behaviour 

of typ© of specimen upto cracking is nearly the same 
regardless of whether it is plain, reinforced longitudina- 
lly or transversely or continuously bound by a helix 
or hoops and longitudinal bars. The specimens behaved 
classically upto cracking and the cracking torque appear 
to depend almost entirely on the geometry of the cross- 
section and the concrete strength and little on the 
amount of reinforcement present. 

Various theories have been proposed to predict 
the ultimate torqjae of R.C.C. members depending upon the 
assumed mode of failure. Rausch, Anderson and co\jan 

based their theories on diagonal tension failure and 

( 21 ) 

assumed elastic analysis to be valid. Lessig considered 
the combined resistance of concrete and reinforcement and 
redistribution of stresses and cleared some of the points 
which otherwise could not be explained by elastic theory. 

main contention was that after first cracking, elastic 
theory is no longer valid. 



Summary of Torsion Theories for 



OTq.sBxd Traes 
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Realizing the importance of reinforced concrete 
members subjected to combined torsion and bending or shear 
or all the three, Wumber of investigatiors ^ have 
studied the interaction, behaviour and failure mechanisms 
R.C.C. members subjected to these loading-^ Number of e5<pre- 
ssions have been advanced to calculate the bending moment 
and ultimate torque. These are based on the equilibirium 
of forces required just prior to failure --echanism. 

Recently, ACT committee on Torsion have advanced 

some tentative recommendations to be incorporated in the 

(26) 

1970 ACI Building code. In a recent paper, HSu & Kemp 
have discussed the background behind such recommendations 
for plain, reinforced and prestressed concrete members 
subjected to pure torsion. The interaction of torsion with 
shear and bending have also been discussed. 

1.3 Scope: The main utility of present work lies in its appl- 
ication ro the practical problems. In the torsion design 
ol machine parts, components of aerospace structures and 
concrete beams, the present analysis assume greater impor- 
tance since materials used for these are generally non- 
linear. However, in order to make fall use of present app- 
roach, further investigations on flanged & other practical 
sections including hollow sections statistical treatment 
of non-homogenity, effect of reinforcement etc. nedd to be 
undertaken. Also, rational theories of inelastic behaviour 
and failures under combined stresses which will help m a 
clear understanding of shear failures, behaviour under- 
combined torsion, bending and shear canaiso be develoned . 
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2, STRESS- STR/iIN BEHAVIOUR: OF STRUCTURAL MATERIALS 

2*1 Introduction 

In this chapter, _ various analytical expressions 
for stress -strain curves are discussed briefly. Out of these 

If m 

few are used as idealisation of stress strain behaviour of 
wide ran^e of structural Materials. Emphasis is laid on 
materials having non-lineat stress strain behaviour. Finally 
by making use of von Mises Octahedral shear stress theory 
and tensile stress strain curves, shear stress strain curves 
are obtained for few structural materials. These curves are 
idealised by Ramberg-Osgood function and the values of the 
parameters are given so that these can be used in predict- 
ing the behaviour of such materials under torque. 

2.2 Preliminaries 

In literature number of expressions have been proposed 
to represent analytically the stress^strain behaviour of 
materials. This approximation of stress- strain curves with 
analytical expressions has received much attention mainly 
because it unifies material behaviour in general by one 
theory, just like Hooke's law in Elasticyty. However, many 
difficulties come ih the way In this sort of representation. 
These arise because We want these expressions to be simple 
enough in order to be of practical use and at the same time 
]^ 0 p 2 resent the given material as accurately as poss- 
ible. The following considerations give a good guide for 


such a selection 
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(i) The analytical expression should be simple and contain 
easily determinable parameters . The number of parameters shouH 
be sufficient to cover’ the stress strain behaviour of wide 
Variety of materials. Two parameters are regarded as the 

most practical. 

(ii) The function should be single valued monotonically 
increasing up to the ultimate stress and continuous. It 
should pass through the origin and hav© a slope there equal 
to the modulas of elasticity. And should give physically 
possible values of stress for all possible values of 
strain or viceversa. 

2 .3 Analytical Expressions 

Polynomial expressions for experimental curves are 
a common practice. Many of the expressions proposed are 
particular cases of the general expressions, given fey 

e=v--+j3^2 + + .... (2.1) 

e + be2 + Ce^ + de^ +.....* . (2-2) 

where e is the strain I ratio of the 

stress, s, to the modulas of elasticity, E. 

( 27) 

Some of the proposed relation with remarks 

about the suitability in use are given m the tabular form 
in Table 'h. 
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2.3.1 Remarks 

Out of the so many analytical expressions proposed, 
only few are suggested as useful. For large strain^ .Prager' s, 
Holmquists and Nadai's & Osgood's expressions are more sui- 
table than Poncelet ' SjRamberg £c Osgood's, Rao & Legget’s 
and Smith & Sidebottom ej^ressions mainly because most of 
these become asymptotic for large strains 4,5 is actually 
observed in practice. First preference IkDweSretr sJobCird'- be 
given to Romberg & Osgood's, Smith & Sidebottoms formulas 
ior small strains & Osgood' s ej^ression for large strain 
Ijcc.iuse of their simplicity & coverage of wide variety of 

nonlinear structural materials. 

2.4 Ramborg-Osgooii .Rdnetion . , . _ 

Ramberg-Osgood^^®^ originally presented the 
functioh in the nondimensional form as 

‘"-=11 + ^ ( -^ )'^ (2.3) 

For the following relation, also known as 
Romberg -Osgood relation, 

fc=^+ (2.4) 

Stress strain constants for various materials 
(29) 


are given in Table 3 
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Table 3 

Typical Stress-strain constant (Approximate) for few 
Materials. 


{ 

FxlO”*^ psi 

B psi 1 

n 

Alloy Steel UTS= 100,100 psi; ! 
Uniform elongation = 2C% ; 

29 

1 

1 

122,400 

25 

Aluminium alloy 2024 -T3 : 
UTS=65,000 psi; uniform 
elongation =15%. 1 

10 

72,300 

10 " 

Aluminium alloy 7075-T6: 

3,000 psi; uniform 
elongation = 1C% 

10 

1 , 

101,200 

1 

10 

( 

Magnesium alloy: UTS=39,000 
i-psi; uniform elongation =15%. 

6.5 

1 

47,500 

10 


in Lh present work, generalisation of Ramberg & 

( ' food's function has been extensively used. For shear 
stress strain in nondimensional form this can be written 

clS 

7 * = ^ 2 . 5 ) 

^here and r are arbitrary constats, 
f, ■«* = i. , Dimensionless shear stress 

< T- 

( 

t 

\ ^ \ , Dimensionless shear strain 

\ h ( may o-ssurne any value recjuired for approximating the 
actu li curve with the above relation. So in a way there 
four parsmeters which provide sufficient freedom to 
choose the values of these paremeters for a good fit # 
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Whenj(_^= this relationship is a ■' "(amberg-Osgood 
function whidh gives a spectrum of curves (frg. 1) ranging 
from elastic to elastoplastic cases for r ranging from 
unity to infinity* For r^l, it is a non‘-linear relation- 
ship. When^= 0/ this relationship corresponds to a para- 
bola of th degree (fig. 2) , By adjusting values of 

experimental shear stress strain curve 
can be approximated using equation- (2.5) 

2.5 Shear Stress Strain Curves 


Till now little has been done to obtain shear 
stro5's strain curve of different structuial materials , alth- 
ough this is very important before -‘ye take up any torsional 
analysis of non-linear materials. This is bocause of the 
experimental difficulties to create such a state pf stress 

and the measurement of the corresponding strain. Buchanan 
obtained the shear stress-strain curves from the experimental 

torque-rotation curve by making use of the following relations. 

' = I? (2.6) 


where T = the applied tor(qae 

r = the radius of the test specimen 
J = the polar moment of inertia of the 
test specimen, 
thc-^ir strain ^ 

0 = twist in the member 
L = length of the monber 


(2.7) 
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The above results are true only for circular cross 
sections. For other sections / these are approximate. 

In the next sub -section, we discuss how to obtain 
shear stress strain curves from unriaxial (tension) stress 
strain diagram by making the ©f the von-Mises Octahedral 
s)iear stress theory . 

2,6 Shear— stress Diagram from Tensile stress strain diagram 
Any non-linear stress strain curwje consists of 
clastic portion Sc plastic portion with the corresponding 
St rains knov/n as elastic and plastic strain respectively. 
Therefore total shear strain 

Wotal ” V Elastic "** Vpiastic 
^ E 



( 2 . 8 ) 


whore G and Gp are modulas of rigidity for elastic and 
plastic portion respectively and are given by 


G = — (2,8.1) 

2(1+Al) 

Ep 

G= (2.8.2) 

2 ( 1 + 

whore subscript P on the constants are for plastic 
portion . 

XX. = *' (29) on the assumption that platic sstraanecaused 
P 

prim inly by shear stresses produces no volume changes. Now 
in order to find a method of determining Ep which is a 
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function of the state of stress, one requirement is *■ that 
the method must give results that agree with the tensile 
stress strain diagram. At any point on the tensile stress 
strain diagram (for plastic strain), the ratio of "Vd 

P 

represents Ep. In order to decide what value of ^should 
be t.aken for finding this ratio, use of the von Mises 
Octahedral sherir stress theory has been made. Since the 
plastic strain are caused by shear stresses , the criteria , 
is based on shear stresses such that Octahedral shear stress 
produced by two different state of stresses x.e. pure shear 
ind )urt' tension, be equal. By this assumption, we obtain 

(2.9) 

where is the effective nomal (tensile) stress 
iul filling the above criteria. 

So in the process of determining Ep, not only the 
value of ^“is obt-j.ined but at the same time a relation 
between shear stress Sc normal tensile stress is also 
obL lined. This relation proves to be of great value in 
gutting the shear stress strain diagram from tensile stress 

strain diagram. 

Substituting (2 .9) , (2 .8 .1) , (2 .8 -2) in (2-8) 

v;i ' gut , 

0.577 r 0.577 v' 

^ total “ G Gp 

_ 0.571-^ + 1.732 C- (2.10) 

G ' 



Thus knowing G and plastic strain & corresponding 
stress for pure tension, the shear stress strain curves 
can he obtained by making use of the aboi^e relation* Shear 
stroc-s-strain diagram have been obtained this 
(fig. 3—6) for few of the structural materials to show the 
practical utility of the above method of conversion. Simil- 
arly/ instc id of Octahedral shear stress, wa could have used 
the P.ijvimum shear stress to be equal for rthe two state 
oC stress mentioned above. But the former ,is preferred. 
2.6.1 Remarks 

The use of such a conversion is open to criticism 
since by doing so we are actually making the maximum 
shear stresses on Octahedral planes due to- two different 
St itc of stress as equal; which othervds^ different 

on other planes. Neveirtheless, in the absence of any shear 
stress strain diagram, this procedure of conversion is 
'worth tloutng. 

Recently a new method of obtaining shearing 
strength of concrete hc».s been developed. The setting of 
Lv sting device and the form of specimen is such that 
th' 'S' avoid other secondary failures and give results 
for pure shear only. This testing device may on further 
niudiiic ition be helpful for obtaining shear stress-strain 
cutVLs of not only concrete but other structural materials. 
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3 . THEORETICAL ANALYSIS 

3.1 I nt roduct i on 

In this chapter the torsion problem for prismatic 
members of unifonn section made up of isotropic and homog- 
enous material has been formulated. The stress strain rel- 
ation IS assumed to be represented by a monotonically 
increasing continuous function *f'. It is diown that in 
general the governing partial differential equation with 
(-he boundary conditions is non-linear in nature and hence 
diificult to solve in closed form. Therefore recourse is 
taken to Rayleigh - Ritz approach of assuming a function 
sat IS lying the boundairy conditions such that the complemen- 
tary energy is minimum. Thus the problem reduces to a mini- 
mization problem which is easier to tackle. 

3.2 Analysis 

Geometry ot Deformation 

Consider a prismatic bar of arbitrary cross-section 
Df lengthLj made up of isotropic, homogeneous, inelastic 
natcriuJ , subjected to a torque T at the ends as shown in 
"iq. 7. The torsion member is assumed to be long enough so 
.hat method of applying the torque at the ends does not 
tli<'’ct the general behaviour at a section away from the 
'lids, . 

The Z axis is taken along the lenght of the bar 



29 


and X, Y axes are taken in the cross-section. The deformat- 
ions of the twisted bar consist of both the rotation and 
warping of the cross-section as assumed by Saint - Venant, 
and may be written as 


u = - Q'zy, 

V = ©zx. 


(3.1) 


© 'p(xyy) 


where n, v and w are the displacements in X, Y and Z 
directions repectively, © is the twist per unit length, 
and y, is the wajg^ing function. 

The non-vanishing strain components are given by 


= © (^-y) 

3X ^ Zz ^ ^x ^ 


= = o C^+x) 

ciY 


(3.2) 


assuming small deformations. By eliminating v/^in equations 

/ 

(3.2), the compatibility equation becomes 


S ■( xz 

Sy 


"Sx 


(3.3) 


Equilibirium Equations 

The equilibirium equation for this problem, in the 
absence of body forces, are 


$ z 


d z 

<3’ (. xz 


y^ = o 


1 XZ yz 


(3 .4) 


O 
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Equations (3.4) express a necessary and sufficient 
condition for the existence of a stress function 0(x,y) 
such that 


Txz ~ 

2 0 
'a y 

= 

fir 


rl 

N 

1! 

ax 

= 


(3.5) 




mm 



where 0 is an unknown function^.'- ’ 

Boundary Conditions 

For the lateral surface of the mevinbei. which is 
free from external fcrces/ the only boundary condition is 

"^xz ^ ^yz 

where 1 and m denote the directions cosines of the angles 
between the outward normal and X and Y axes respectvely. 

1 ss dy/ds / m = - dx/ds 

where s is the avc parameter. 

Substituting (3.5) in (3.6), we get 

ll SI ^ ^ ^ ,0 ( 3 . 7 ) 

dy ds "Sx ds ds 

This is satisfied if 0 is constant on the boundary. 
Since the stresses are given by the partial derivatives 
of 0 , the boundary condition may be written as 

0 = 0 on the boundary (3.8) 

i 

At the ends of the bar due to the assumed stress 
distribution, resultant forces in X and Y direction are 
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zero . However the couple foimed by these forces is equal to 
the applied torque there and is given by 


T 


j[ - y axdy 

Area 

('ll X dx ay - Y dx ay 

A A 

2JJ dx dy 


A 


( 3 . 9 ) 


Constitutive Equations 

The relation between the resultant shear-stress X 
and resultant shear strain Y can be written in the follow- 
ing general non-dirnensional form 


V 

Y* 


= f (^*) 


Vo 



( 3 . 10 ) 


where T== yTxz ‘’'Xyz 

and T^rnay assume any value required for curve fitting 

and non-dimensionlization . The function '£ is a single 
valued monotonically increasing continuous function, which 
may be linear or non-linear. 

In order to express the inelastic stress strain 
relationship, we make use of Hencky's total strain theory- 
on the following assumptions, 

1 

(I) Wihciple axes of stress and strairi coincide. 

(II) Mohr’s cirdle diagrams of stress ahd strain are similar 
at any stage m the inelastic deformation. 
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we get the following relation, 

■3xa = Xxz 

\xz Yyz 


(3.11) 


Substituting equation (3.11) in (3.10), we obtain 


Vxz = Vo To 0y f (T*) r * 

= -VoXo^j^ f ('^) * 


( 3 . 12 ) 


Equation (3.12) when substituted in equation (3.3), 
leads to a non-linear partial differential equation for the 
stress function 0 , with the boundary condition that 0 
should be zero along the boundary. However, this equation 
becomes too complex for a general function ' f ' . Hence 
recourse is taken to the Rayleigh Ritz method'* 

3,3 Rayleigh Ritz Method 

By making use of the principle of minimum compleme- 
ntary energy which holds good even for non-linear materials, 
this method becomes a powerful tool for solving boundary 
value problem in the field of elasticity^ especially when 
the governing differential equations arr difficult to solve. 
Here a function in the series form with some unknown coe- 
fficients (i = 1,2, ...., n) is assumed 5uch that boun- 

dary conditions are satisfied. These unknown coefficients 
are determined by minimizing the complement airy potential 
with respect to each A^, leading to n equations in n 
unknowns. By solving these, 0 is completely known and 



33 


thus for a given value of $2? , the torque ahd the stresses 
can be calculated using equations (3 ,9) and (3,5) . 

3,3.1 Completnentary Potential 

The complementary potential ~JT is given by 
TT = (3.13) 

where = Total complementary energy in the number 

= J d (Vol.) (3.14) 

Vol. 

where = Complementary energy per unit volume, 

= j '/oTo ^ (3,15) 

and Wg = the total complementary work done by the 
external torque acting through the angle 
of twist 9. 

= 2 L J j" © 52^ dx dy (3 .16) 

A 

Substituting equations (3.14), (3.15) and (3.16) in (3,13), 
we get, 

Tf = L Vo-To j j{ ^T* ) <1 T* - ax ay O . n) 

A 
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4 . APPLICATIONS 


4*1 Introduction 

In order to demonstrate the applicability of the 
theoretical analysis discussed in previous chapter, two 
casfes> rectangular and circular cross-sections are consi- 
dered in this chapter. Smilar procedure can be extended 
to other practical cross-sections. From stress strain 
relation stand point, only generalised Rarriberg-Osgood 
function is replaced for the general function ‘f*. This 
IS done with the view that such a function with suffic- 
ient member of parameters covers almost all non-linear 
materials and hence the present analysis can be applied 
to a wide variety of materials . Although here also, it 
can be extended for some other stress strain function ’ f ’ 
which may be more suited to any particular material. 

4.2 Idealised Non-linear Stress Strain Law and Comple- 
mentary Potential 

Generalised Ramberg-Osgood function can be applied 
to a geheral class of non-linear materials and is given by 

y* = 

For this relationship the Complementary potential Tf 


can be written as 

-iT = I 0 

A 






2 

+ 




m 


dx dy (4.2) 
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where m = — 2 — 

Another curve of interest is of the fom 
Y* = Sinh X* (4.3) 

Which has been extensively studied by Smith and 
Sidebottom (ll) for the study of circular, rectangular, 
and triangular cross-section with particular te:;erence 
to steel members . This curve has been appoximated by the 
Ramberg - Osgood function by taking the first two terms 
and the results have been discussed in section 7« 

4.3 Circular Cross-Section 

For a circular cross-section of radius 'a' and 
length L, the function 0 can be assumed in general as 
follows, 

P q , 

^ (x,y) = ( 21 H ( x2+y2_a2) 

1=0 j =0 

We retain only the terms with coefficient ^ 20 '^ o 2 - 

Because of symmetry A 2 Q = Aq 2 

^ (x,y) = [a.qo"*''^ 20 (x^+y^)\ (x^+y^-a^') (4.4) 

to 

or in polar co-ordinates 

0 ; (r, ©) = (A+Br2) (r2-a2) s (4.5) 

to 

where A = Aqq 
B = A 20 

(^)+C^)^= 4 r2 


(4.6) 



where R 
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/ 


r 

3. 

Aj^ = As 
» B 

Substituting (4.5) ^ (4.6) m (4.2) and integrating 
we obtain 



where h = 2a 

To minimize this ccmplementary pbtential in-order 
to determine the value of the unknown constants A and B, 
we equate the partial derivatives of the complementary 
potential w-.r.t. A and B each equal tc reyo 


3lT 

A 


0 



4m f j^2m+l^^^_l_2BiR^-.B^) 


2m- 1 


dR 


©h 

+ h — ^ O (4.8) 

^ o 

and *^J!L =s 0 

dB 

or ^ ^ ^ (2R^-1) dR 

I 0h 

12 Vo 


4- 


(4.9) 



37 


Tvhese two simultaneous nonlinear integral equations 
are to be solved numerically for and by the method 
described in the next chapter. 

For a given value of 0, the torque T is given by 

To 2a “ 

and maximum is given by 

max ~ ^ (a^+B^) (4.11) 

4.4 Rectangular Cross-section 

For the case of rectangular Cross-section, of 
width 2a, depth 2b and length L, stress function 0 can 
be assumed in general form as 


^(x,y) 


p q 

^ z: 

1=0 j =0 


x^y 




(x^-a^) (y^-b^) 

(4.12.1) 


q 


or 


Ltrx 


A Cos Cos ^ (4.12.2 1 

i-i,3. 2a 2b 


Retaining only the terms with coefficient A^^, A^2 and 
A 2 q in equation (4.12.1), we get 


^ = (Aoo+-^ 20 x"+A^-y2) (x2-a2) ( y 2 _ b 2 ) 


(A+Bx2+Cy2) (x2_a2) (y2_b2) 


( 4 . 13 ) 
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(^) ‘= 4^2 [[k^X^ (Y^-l) ^ (2 x2_i) +K^CjY^ ~ ^ 

+ y2(x2_1)2 ^A^+B^x 2^^K^(2Y^-1)12 "j 

(4.14) 


where K = b/a 

X = x/a 

Y = y/b 

Substituting (4.13) and (4.14) m (4,2), we get 

2.84444(K^+1) K^A^+(0 .893968 

+ 0.135450) Bi +(0.893968+0.812699 K^) 
k7cJ + (1.13778 K^+ 0.812699) K^A^b^ 

+ (1.13778 + 0.812699 K^) k5a3_C^+0.162540 
(k2+ 1 ) k5b^c^-+-^ J'^|(||)^+(^)^J ’^dx dy 
•' O o 

- % (SA^ + Bj^+ k2 c^)J 

( 4 . 15 ) 

where A^ = Aa^ 

= Ba^ 

= Ca^ 
h = 2a 

To determine A, B and C, we equate the partial 
derivatives of (4.14) w.r.t. A, B and C equal to zero. 

Thus we get 

dir 

d> A 






0 
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and 


or ^£5.68888 (k2+i) (1.13778 k2+0.812699)B, 

+ k2(i. 13778 + 0.812699 K^) Ci"! 

+ 4 '"-! K 2 -»- 2 jJg 2 x 2 (Y 2 -i,^^^';;Jg^ 

4- ^ +y2{x2-1)2£ii^h.BjX2+Ci(2y2-1)}2J”’-' 

[k2x2 (y2.i) 2 £ai«^(2x2.i, +K2c^y2 7 
+ Y^(X^-1)2 £aj8Bjx2+k2c7(2y 2-1)JJ « dY 

' 9Y^ ~ ° (4.16) 

btf 

= 0 

o B 

^£(1.1787936 K2+0.270900) 8 , 3 .( 1.13773 k2+ 

+ 0*812699) + k2^0. 162540 (k2+i) 

, 4m-i 7^2in-2 f^V 

J|[k2x2(y2_1)2£7,^^B,{2x2_i, 

+K^C^Y^ } ^ + y 2 (x 2 _i )2 +B 2 ,j, 2 p 

■« ClY 2 j ( 2 x 2 - 1 ) + y 2 (x 2 -i) 2 £ 3 ,^ 3 , 3^,.2 
+k2c, (2y2_i) J x 2~[ dX dY 



o 


(4.17) 


3^77 

3 ^ ~ ‘ 2 > 
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or 


^ [^(1.1787936+0.270900 K^) k2c^+ (1.13778 
+ 0.812699 k2) A;j_ + 0*162540 (K^+l) B J 

^ ^m-l j^2m-2|'£^^2x2(Y2_l)2^Ai+B3_(2x2-l) 

+ k2 C]^y2 ^ 2" % y2(x2-1)2 Ai+B3_x2+k2c^ 

(2X2-,!) ^ 

+ k2Ci y2 *7 y 2 + Y2(x2-l)^r A2^+B3_x2 +k2C3_ 


(y2-1) 2x2-1) 



dX dY 


©h 

45^0 


(4.18) 

7*gain these three non-linear equations (4.16)^ 
(4.17) and (4.18) are to be solved numerically for A 2 _,B^ 
and by the method given in the next chapter. 

For a given 9 , torque T is given by 


Xo 4ab 2a 


45 


(7.3_+B^4C^) 


(4.19) 


and maximum 'X' given as 


X 


★ 

max 


2 (A^+B^) 


(4.20) 
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5. NUMERICAL TECHNIQUES USED 

5,1 Introduction 

Problem encountered in the previous chapter can be 
attempted either analytically or numerically. Analytical 
solution/ however, does not seem possible here, firstly, 
because the simultaneous equations are of non~linear in 
nature and secondly the integral involved in each of the 
equations can not be evaluated for non-integral values of 
'm' in general, unless the values of A,B, C are Tcnown. There- 
fore a numerical technique is suggested and used for the 
solution of the abo^re problem, A particular numerical proc- 
edure is the usual Newton' s Iterative technique to solve 
the system of non-linear equations. At every iteration first 
the integral involved in each of the equations is evaluated 
by the Simprous Rule, for the values of A,B,C obtained at 
previous iteration. To start with the method, suitable values 
of A,B,C are chosen. Once the integrals are evaluated, the 
system of nonlinear integral equations reduces to system 
of non-linear equations and therefore can be solved by 
Newton' s method. 

5.2 Simpson's Rule for Numerical Integration 

For definite single integration, the simpson's rule 

x_4-nh 

^ It. 

Jydx = - (yQ+4y2+2y2+4y2+2y4+4yg+ ) (5.1) 



• mm* 


denote the ordinates at 


where y^ 

Xq,x, i.,... etc and h is interval length. 

For double integration, the extension of the above 
rule known as Extended Simpson's Rule is as follows, 

X +2h y +2K 

J I dx = — [z^^+Zo2+Z22+22o+4 

>^o "S. 

(Z01+Z12+Z21Z10) 

(5,2) 

where denote the functional values 

at grid points and h, k are grid lengths 



g 

Tie coeff. of Z are shown above in the figure 
By adding any number of such unit blocks the whole domain 
of integration can be covered. 

5,3 Solution of Simultaneous Non-linear Equations 
(Nowton' s Method) 


Let the three equations be 
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f(x,y, 2 ) = 0 

g(x,y,z) = O (5.3) 

h(x,y/z) = O 

Let (xQ^yQ/ZQ) be an approximation to the solution 
( ) of equation (5,3) . 

Assuming that f,g and h are differentiable, we 
expand f(x,y, z), g(x,y,z), h(x,y, z) abottb 

by Taylor' s Series for functions of three variables as 
f(x,y,z> = £(xo,y^,ZQ) + 

+ fy(xo,y^^,ZQ) (y-yQ)+ 

(z-Zq) + 

g(x,y,z) = g(xo,y^,z^)+g^(x^,yp^z^) (x-Xq) 

+ gy(xQ,yo,Zo) (y-yo)+ 

( z-Zq) + 

( 5 . 4 ) 

h(x/y,z) = h(xo/yo/Zo^ ^x^^o'^o'^o^ 

+ hy(xo,y^,Zo) (y-yo)+ h^{x^,Yo>^o'^ 

+ 

Further, assuming that sufficiently close 

to the solution ( / 't ) r higher order terms can be 
neglected. Therefore equating the linear terms to zero. 


we obtain 
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(x-Xq) 4gy(Y-yQ) 

(x-Xq) +h^(y-yo) +h2(z-ZQ) 


f 

g (5,5) 
h 


where it is understood that all functions and their partial 
derivative are evaluated at (Xq/Yq/Zq) , We then expect that 
the new solution (x 2 /y]_/Z 2 ) of (5.5) will be closer to the 
solution ( / '^ ) than (xq/Y^/Zq) , Then by Cramer's 

Tjrule 


xi-x^ 


- f 

- f 

- h 


gy 

hy h^ 


1- f 

- g 

- h 



f 

Z 


5z 


h 


z 


"x 


g 


X 


h 



Yl-Yo = 


^x 

-f 

fz 

Px 

-g 

gz 

^x 

-h 



J(f/g/h) 


Z-Zq = 


fx fy -f 

■gx % -g 

hx hy -h 


J Cf,g,h) 


J(f,g,h) 
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Provided J(f,g,h) at ^ 

The solution (x 2 _ 7 yj_/ Z 2 _) of this system now provides a 
new approximation to ( 2r r Y) / ) 

In this way we generate successive approximation till 
the convergence is achieved^ 


Xx 


■i+l=^x + 


'i+l 


= yi 


'i+l 


Zi + 


- f 

fy 

f 

z 

- g 

fy 

gz 

- h 

hy 


J(f,g 

/h) 


^x 

-f 


5x 

-g 


h 

z 

-h 

h 

z 

J(f 

/g/h) 


■' fx 

fy 

-f 

^x 

% 

-g 

: h^ 


-h 


J(f.g.h) 


5.3,1 Convergence of Newton's Method 

The above method when converges, it converges 

( 1 \ 

quadrat! cal ly. It has been shown that this method 

will converge under the following sufficient conditions? 
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( 1) £,g/h and all thexr derivative through second order 
are continuous and bounded in a region R containing 

r; ) 

(2) The Jacobian J (f,g/h) does not vanish In R, 

(3) The initial approximation is sufficiently close 

to the root ^ ^ 

5.4 Selection of Sub-inteirvals and Convergence Limit 

For stable and sufficiently accurate results / it 
IS necessary to select suitable no. of sub-intervals for 
the numerical integration. Although efficient computation 
demands the no. of intervals to be as small as possible/ 
it IS rather difficult to put analytically any limit on 
it to ensure stable and accurate results . These are best 
selected by meams of trials. For the present problem/ 
twenty subintervals have been found to be adequate* 

Convergence limit which control the no. of itera- 
tions are selected depending upon how the results are 
going i to be used. For plotting purposes twenty iterations 
or convergence limit 10*"^ whichever is readhed earlier 
is quite sufficient. 
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6. ULTIMATE STRENGTH OF PLAIN CONCRETE IN TORSION 

6.1 Introduction 

For the design of concrete structures prediction 
ultimate strength is as important as the behaviour - frc~ 
begining to failure . In this section a methodology^’ has 
been presente for obtaining the ultimate strength of plain 
concrete in torsion from uniaxial test results. Sinh cuarve 
has been used to represent the shear stress strain relat- 
ion. Result for rectangular cross section have been 
obtained and compared with other theories and experimental 
results. Such a methodology can be extended to other sections 
and stress-strain, functions. 

6.2 Tensile Strength of Concrete 

The shape of the stress strain curve in pure sheaY 

( 71 — V) IS similar to the shape of the stress-strain curve 

( 29) 

in uniaxial tension ( o- - (L ) . The conversion factors 

depend upon the failure theory applicable fort the material. 

For example, the octahedral shear-stress theory which is 

(29) 

applicable to steel, gives 

( — and ^ (6.1) 

O' 

For concrete, the above expressions are very conservative 

Mohr's failure theory is applicable for converting the 

maximum shear stress to equivalent uniaxial tensile failure 

stress. This can be written as 

X = f^ where f^ = ultimate tensile stress of 

the material 
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The value of ixfusing a straight lime simplification 


(13) 


of Mohr's theory is 


f 

c 


, while Iyengar 


(33) 


c ^ -"t 2fl - f' 

recommend an approximate value of — S v/here f = 

2f^ ^ 

c 

ultimate cylindrical strength of concrete. McHen3ry and 
(34) 

Kami also give values of from their experimental tests 

ovaries from 0.85 to 0.95. The value of f^ for design 

(13) 


purposes given by Hsu 


for stone concretes is 


f 

t 


= p - '/ 1 


where p' = 5.0 


( 6 . 2 ) 


However, in analysing test results ^may vary between 4.5 to 

(35) 

6.5 . Tor high strength concretes Navaratnarajah uses 


the following relationship 


f = 0.68 (f‘) 
t c 


'S 


(6.3) 


The tensile strain for a given value of shearing strss 

, (33) 

can be written as 

V = 2(l+/JL)e (6.4) 

where IX. = poisson' s ratio, which is assximed to be 

constant over the whole range of loadin j, with a value 
(13) 

of 0.17 . The normal stress-strain relation («?" - 6 ) 

from (X - V" ) relationship is therefore given for the Sinh 
relationship (equation 4,3) as 


1 + ^ 4 ) e <r 

— = Sinh 


V 
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At ultimate, the value of cr'is ^ with a value of 

the value oft;is^ 4 The value ^f^from direct tension 

stress-strain curves is reported in the litesrature to be 

—6 

around 150 x 10 by Hsu, 200 -* 400 x 10 by Evans and 
Marathe''^^^ , 400 x 10 “^ by Mirza^^”^^ « Further research is 
needs to fi,nd the parameters influencing this value, 

6.3 Ultimate Strength Methodology 

With the above considerations a methodology for 
finding the ultimate strength, angle of twist per unit 
length O of prismatic plain concrete members_in pure tor- 
sion, IS herein presented using the Sinh curve idealisation 
for behaviour of plain concrete* This is applied to the 
case of rectangular beams . Now 


V 


Sinh _ 
TI* 


where are to be determined in terms of f ‘ , f * 

fy ^ c t 


the origin, where 


It can be readily seen that 


c 

= G, 


at 


V 


' (1 


Gq = 2 ' ("f = shear modulus 


( 6 . 6 ) 


Equation 4,3 can be approximated as 
- _ X 1 

by neglecting higher order terms in the Sinh expansion 
Using the above ej<pressions, the following relation is 


obtained 
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To 






The value of E in tension is taken as the same as in com- 

(35) (->3) 

pression after Todd; Oladapo and Hsu'' 

The value of E can be written as =57/500 /f} 

f ^ C 


for stone concretes. Using t = 250 x 10' 


^ = 5 the value of can be written as c.^,= 
Example 


-6 

_ T 


= 0.85 


3.8 


Using the results of Smith and Sidebottom ^^*^^ , for 


0/a of 1.0 the value of 
Q 2a /~7 

by 


54 and the ultimate torque ^ is given 


T 


8a^b 


= K.f 


where = 0.254. 

The coefficient 0.254 in the above expression is only a 
probable value and may vary between 0.23 to 0.275. For 
clastic analysis neglecting the presence of oompression 
stress the coefficient is 0.208 while it is 0,33 by rigid 
plastic theory. 
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%t, RESULTS 

7 . 1 Summary 

Firstly a stress Sanction satisfying the boundary 
conditions was assumed. It dohtained two and three unknown 
parameters for rectangular and circular cross-section 
respectively. Substituion of stress function in the comple- 
mentary potential derived for stress strain function given 
by Ramberg & Osgood, resulred finally in making the 
complementary potential a function of the unknown parame- 
ters. These were determined by making use of the principle 
of minimum complementary potential. This minimisation w.r.t, 
each of the unknown parameters led to a system of non— linear 
integral equation, which were solved numerically on IBM 7044 
for different values of twist by the numerical technique 
described in section S. The computer program is given in 
the appendix. 

Once the stress function is known, T* Vs 9* and 
max ® curves are obtained making use of equations 
4,10, 4.11 for circular cross-Section and equations 4.19, 
4.20 for rectangular cross-section. All this done on the 
assumption that the values of the unknown parameters of 
the idealised shear stress strain curve are known. 

1^.2 Circular Cross-section 

Results are obtained for different materials. For 
concrete the value of r = 1-10 to 1.16 for S =■ 1 was 
suggested by Marshall & Tembe, while the parabolic curve 
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(r = 2, = O) was used by Miyambo, in their analysis of 

circular cross-section in pure forsion. Curve for T* Vs 9* 
fear different: values of r and ^ for concrete are shown 
in fig. 8. 

For other non-linear materials like Alluminium, 
Magnesium and Steel, the T* Vs 9* curve is shown in fig. 9. 
The values of parameter of Ramberg & Osgood function for 
these materials are given in section 2 . 

7.2.1 Comparison 

Results are compared with Buchana^t’s experimental 

f 

results for beams Cj^, C 2 and C^ . The 'C - x curves obtained 
from experimental results are approximated by hyperbolic Sinh 
function and Ramberg-Osgood function. The comparison for 
T Vs 9 and T* Vs for these functions is done in fig. 10 
and fig 11 respectively. 

Besides this, results are also compared with Smith 
& Sidebottom. The hyperbolic Sinh stress strain relation used 
by them is approximated by the Ramberg Osgoods function by 
taking the first two term of the expansion as given below 



r 


3 i .e . m = 2 
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For these values of V and £ , results for T* Vs ©* 

and <- * Vs ©* are compared in fig, 12. Results for T*, 
max 

X* for various values of r and“^ are given in Tahulalar 

form in Tables 4 and 5. 

6.3 Rectangular Cross-»section 

For this case curves for T* Vs 9 for different 
values of r^ = 0 for concrete are presented for b/a = 1.0 
1*5, 2.0 in figures 13,14 and 15 respectively. 

T* Vs 9* curves for aluminium, magnesium Steel and 
reinforced plastics are shown in fig. 16. 

6.3.1 Comparison. 

For rectangular case also, results are ccsmpared 
with Buchanan's e5<perimental results. The T'- rY curve was 
idealized by Sinh & Ramberg Osgood functions. The curves *■ 
for T* Vs ©* are compared in fig. 17. Ramberg-Osgood func- 
tion apt)roximation give better results than Sinh function. 
Comparisons with Sidebottom & Smith' tidal 

K 

results are shown in figs. 18 & 19, for T* Vs 9* and f 
Vs 9* respectively. The Sinh function is approximated by 
taking first two terms of the exparison as given earlier. 

Results for T* , -f' * for various values of 9"^ 
r, ^ and b/a are given in tabular form in Tables 6-9. 
Ultimate Strength 

The values of by the proposed method (section 6) 
are shown in Table 10 and compared with values obtained by 
elastic and plastic analysis. From Table 10, it can be seen 
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that the present analysis gives a value of t ultimate streng- 

th/ about IQYo of that by plastic analysis. This is in comfo- 

i 3 8 ) 

rmity to the findings of Colling who' proposed the Nadai's 

sand heap analogy with the value of the stress to be 3.5 
To check the validity, of above ultimate torque relations 
with e3<perimental result of rectangular beams, compabison is 
made in Table 'll. The value of f'^ is taken to be the value 
obtained by direct tests (where available) , on the maximum 
of the value obtained from the expressions given in this 

^ 38 ) 

sedtion. In die tests reported by Collins, the value of 
f^ IS taken as 0,56 times the value from split cylinder tests. 

From Table 11 it is seen that the ultimate torque 
obtained agrees very closely, when the actual value of 
is known and is conservative when p, = 5.0 is used. The value 

I 

of twist at ultimate is not compared, as this is very much 
dependent on the value of ^ assumed in the general esq^ressions . 
7 .4 Discussion 

Results have been presented for five structural mater- 
ials. However, the comparison with the experimental results 
is done for concrete only since these were the only avail- 
able. For both cases, circular and rectangular cross-sections, 
two curves have been used for approximating the given X ~ Y 
curve. In both cases, Ramberg -Osgood function give better 
results than Sinh function, although both these slightly 
over estimate the torque. More experimental results for 
different materials need to be obtained and compared before 
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a definite conclusion can be made as to which function give 
better results for a particular material. 

Comparison with the analytical results of Smith & 
Sidebottom, obtained after approximating the Sinh function 
by Ramberg Osgood function, show good agreement particularly 
for circular case. A point to be worth noting is that ana- 
lytical results are exact for circular case as given by 
Smith and Sidebottom. For rectangular cross— section at one 

of the places in derivation Cosh -c is replac^djoy + 

4 ^ 

n^ , m = 0.47185 & n = 0.057442. This enstt:^?^ maximum f|rror 

of 1.5% if 0 . 

Proposed ultimate strength methodology for plain 
concrete on comparison with clastic & plastic theories and 
experimental results show good agreement. 

However, more research is needed on the tensile 
stress -strain relationship for predicting the ultimate twisty 
as this IS very much dependent on the value of ultimate 


strain 
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Table 4. T"i' for Various Values of 0* for 
Circular Cross-section 
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9 . SUMMARY AND CONCLUSIONS 

In - tla±^ 3di's sertation, a method based on complemen- 
tary energy have been presented for obtaining the torque 
Vs rotation and shear stress Vs rotation curves for pris- 
matic members - The method can be applied to prismatic 
members of any arbitrary cross-section with isotropic, 
homogeneous material properties, provided shear stress 
strain curve of the material is known. 

Shear stress strain curves for structural materials 
at present are obtained either from experimental torque Vs 
rotation curves for circular beams or can be obtained 
from tesnile stress strain curves of the material. But 
none of these methodsgive curves which represent actual 
shear stress strain curve in pure uniform shear. More 
research need to be undertaken for obtaining these curves 
in pure shear directly to replace the preaert methods . For 
this, modifications in the recentlyy developed technique 
of obtaining shearing strength of concrete in pure shear 
can also be considered. 

Idealisation of shear stress strain curves is 
done by generalised Ramberg Osgood, hyt>erbolic Sinh and 
polynomial functions. Ramberg Osgood function is found 
to cover the widest range of material behaviour and hence 
used exclusively. However, there is a limit beyond which 
curve fitting can not be improved. Further attempts lead 
to spoilling of previous fitting in either the initial 



stages or in the final stages of the curve. The possibility 
of using two curves can be thought of as the possible 
answer. The other alternative is to findout some other 
function more suited to any particular material and then 
extend the present technique for that function. 

Results for prismatic members of circular and 
rectangular cross-sections have been presented for mater- 
ials like concrete, aluminium, steel magnisium and reinfor- 
ced plastics. The number of terms in the stress functions 
have been kept to two or three only. More no. of terms 
complicates the analysis and are not practicable. Membra- 
ne and sand tieap analogies suggest the possibility of 
making stress function to depend upon material properties 
as well besides satisfying the boundary conditions. 

Comparison of results with the limited ejq^erieme- 
ntal and analytical results show good agreement. It is 
observed that the Sinh curve for concrete gives value 
of ultimate torque which are in agreement with experimen- 
tal results for reccangular sections. More knowledge on 
the stress-strain curves in uniaxial tension for differ- 
ent concrete is needed for getting the value of twist at 


at ultimate. 
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APPENDIX 

DIGITAL COMPUTER PROGRAM 

A.l Introduction. 

A digital computer program is ffoittan to 
handle nurherical procedure outlined in section 5. The program 
consists of main program and one subroutine. It is quite 
general, and can be used to obtain dimensionless torque 
and dimensionless shear stress corresponding to any value 
of dimensionless rotation. The program covers almost all 
structural materials provided suitable values of 'r' and 
’ ' are given. For the understanding of overall logic and 

linkage of the program, the flow chart is given. Program 
for rectangular cross-section is given at the end. For 
circular cross-section, a similar program will work. The 
terms for input and output variables used in the program 
are given below. Other symbolic names used in the program 
are self e2q3lanat03ry . 

A. 2 Definition of Variables 
(a) Input Variables 

A,B,C = Values of the constants to start the 
Iteration . 

b/a ratio for rectangular section 

Value of m ( = ) 

2 

Number of divisions for numerical 
integration 
Value of ' * 


RAT (I) ,RATIO= 

EMM(J),EM 

ND 

DELTA 
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(b) Output Variables 

Variables which are same as input variables 

are not defined here, the remaining variables 
are as follows 

I IM = Number of Itration 
THETA ^Dimensionless twist 
To r qu e= Dim en s 1 on 1 e s s To r qu e 

T0U(Max)~ Dimensionle ss Shrear Stress (maxinium) 






5 . 


i-' 






'|fi?i 






^ r) ^oo non 


T0'-3Ic^' OF ■^ECT/if'GJL-''^ CPOSS-SECT IC^' DF ^^^'-LIMFAo ‘iATEn>lAL 

3 IVEAGIO’'; U ( 3C ) ,PAT ( F ) ,E'i| (6) 

FEAOIG > ( '^AT ( i ' , 1 = 1 ,3 ; , ( E'^‘/( J ) , J=l .6 ) jf'O 
•10 rOR||-~t] >'«5rfa,?,2x»6F5.2»^X?I3) 

READll , \ ,B,C sDEL” 
i: -O'i'i^T ( 1/ j FIO.A 1 

P I N T 1 2 

12 F0PfAT( ///40X,5 oH TORSI j 1 OF FECT A^s JLAF 3AR OF 
iNCriLlfJEAF I.ATEPIAL/) 

D32^I 1=1 ,5 
D02AJ1=1 ,6 
i^ATIO = RAT( li) 

F l = RATlO-^-RATIO 
E,M = E.-i '( Jl) 

^r<I,vTl3,'^ATlO,Ci'',, delta 

12 for v'rtT ( 37 < ,7H d/^ = , F6 . 3 , 15/ ,4HEM= » F5 . 2 , 1 5X , 7HDELTA = ,13) 
EM il = E’'-l. 

GE' RATION OF CcEFFICIEvT F'^P CIMpso^jS InTCGpATIQN FopMUlA 

f;i=i4D-rl 

DCx'^' ix=2 ,ra ,2 

o'( IX )=A,on 

IF{ IX.EO.I'jD) go to 14 
U( IX+1 )=2,00 

14 CONTINUE 
U(1)=1,00 
b {N1 ) =1 , 00 
D IV=ND 

SEGMT = 1 . 00/DI V 
PRINTI'"’ 

15 Format (/ 3oX,5HTHETA,13X,4H I TRN,15X,1HA,1?X,1HB,19X,1HC/) 
D023L=2 ,1001,5 

l) = L-l 
S = N 

D019i«1=l ,20 

INTEGRAL EVALUATION 

X=0.00 

V = 0.0 

VV=0.0 

VVV=0.0 

V11=0.0 

V12=0.0 

V21=0.0 



\/22 = 0.0 
V31=0.0 
V32=-0.0 
v'Vll^O.O 
VVl2=o.O 
VV21=0.0 
VV22-0.0 
V\/31 = 0.0 
\/\/32 = 0.0 
VVll = 0.0 
WV 12 = 0*0 

vV\/21 = o,o 
\,'VV22= 0.0 
v^V31=o.O 
vV'V32 = 0.0 
roi7i = i,fii 
X2 = X 'X 
X3=/2-l, 

X‘+ = X3*X3 
X5=2. XO”!. 

!■' = 0 • 0 
0 = 0 • C 
R = 0.0 
FI 1 = 0.0 
P12 = 0.0 
F2i=0.0 
F22 = 0.0 
F31 = U.0 
?32 = 0.0 
111=0.0 
012 = 0.0 
u '21 = 0.0 

u22 = 0.0 
031=0.0 
032 = 0.0 
R1 1=0.0 
P12 = 0.0 
i32 1 = 0.0 
Pv2 2 = 0. 0 
P31 = 0.0 
R3 2 = 0 . 0 
Y = 0.0 

D016J=1,N1 

Y2=y*Y 

Y3=y2-1. 

Y4 = Y3*''3 
Y5 = 2.-:^Y2-1. 

CALLSUBI ( XjV »A,b,C,X2,X3,X4,X5,Y2,r3,YA,v'5,Efl,Rl,F7,F8,F9,Fl0» 

1 Fll»Fl2,Fl3,Fl4,Fl5,Fi6,Fl7,Fl8,Flo,F2C,F21,F22,F23,F24,F25,F2 
2F27 ) 

P = P + F7v'J( j ) 

0 = G + F6-'!-'J ( J ) 

R = R + F9"--U( J) 

Pll=Pll+rlO*U£J) 

P12 = P12 + F11:^0 ( J ) 



P21 = P21+Fl2-'tU( J) 
F22=P22+F13*U ( J ) 
P31=P31+F14*U ( J ) 
P32=^^32+F15»U ( J ) 

Q1 l-Qll+F { J ) 
Q12=Q12 + F17 -!:u ( J ) 
Q21=Q<:l + FlS^t-U(J) 
2)22=Q22 + F19-:^U { J ) 
C!31=S31-rF20”-U( J) 

Q32 = Q32 + F21*u{ J ) 
Pll=Rll + F22-'tU( J ) 
R12=R12 + F23"J ( J ) 
P21=R2l4-F24*U ( J ) 

R22 = R22+F25-"-U ( J ) 

R31 = R31 + f26;:-U ( J ) 

R32 = R32-hF27-;!-U ( J ) 

Y = Y + SEG.''1T 

16 CONilNUZ 

V = v+P-;( j( I ) 

VV = Vv'+Q!^U ( I ) 

VVV = VVV4-R*'j{ I ) 
vll = V] 1 + P11^'J( I ) 

V12 = V124-Di2ifU{ I ) 

V21=V21 + P21^U ( I ) 
722=^V22+P22'>% { I ) 

V31 = V3H-P31vu ( I ) 

' y3: = V32+P32-''-U( I ) 

Vyil-YVll+Qll-ifjt I ) 
Vyl2 = VVl2-Ql2-''U ( I ) 
v'V21 = VV2] +G21-:--U ( I ) 
Vy22 = VV2<^+G22-:cU ( I ) 
/V31 = VV31+Q31rcU ( j ) 

n/'v 32 = VV32+Q32-'«-U ( I ) 
VVV11 = VVV11 + R11'»U( I ) 
Vyvl2 = VVVl2 + R12 !.U( I ) 

yvy21 = vV’'21+R21*U( I ) 
yVV22 = vVV22+R22-:!U{ I ) 
VVV31 = VV\/31 + R31 xU( I ) 
Vyy32 = VVV32+R32^:u( I ) 
X = X-^SEG'''T 

17 C3NTINUE 
DD=GEGMT^3EGMT/9. 
y= y-^'GD 

yv=\, 7*00 
V\/V = VV 7* DO 
71 1 = V11*DD 
V12 = V12'-DD 
721=721*00 
V22=v22*DD 
731=731*00 
732=732*00 
7V11 = 7''11*DD 
V7l2 = V7l2-xDD 
V721=V721*DD 
y7 2 2 = V72 2 *DD 



V\/31 = VV31-;<-DD 
\/^3«i = VV32*DD 
yvvii=vv'ai!:DD 
'yVVl2 = VVV 12^DD 
VVV21 = VVV21 !v-DD 
VVV^2=VVV22*DD 
yVV31 = v'VA/31 !^DD 
'v yV:i2 = VVV3^ -DD 


C 

c 50 Lcric''J CF mqnlpieap Equations ey riE-.Tcis method 
c 

TE,'ll = 4.-^-*Ef''ll 
TF'i2 = 4.-K^Ef'i 
Tl = EK-11/2 . 

TR12 = P1-^^. (EM-l. ) 


2 EVALJATIO\ OF FuMIOmAL VALUES A.iiD THEIP P^PTI^L DERIVATIVES 

D=DELTA 

/-F = ( ( 5.6S.363-'> ( Rl + 1. } 'f-H+ ( 1 . 1 3773 :'-p 1 + 0 . 3 1 2o R = ) 
l + Rl-"- ( 1.1 3 7'’3+0» ?1Z6RR“P1 ) - C ) -'-D/32 . ^ + TE.‘'l ;-TRl2 •< V-S/o . 

AG= ( ( ( 1 .TSTSSb / pL-* V. 27'J''CC ) f'5+( 1. l3''78"Pl-^0.? 126RR ) -“-A 
1 -El ' (v. 162340'' C?i + :. )-'-C) ' '^D/32. ) +TEU 1 '-TP 1 2 “ VV- S/ 45 . 

Hi:=( ( (1.7E7925 + v.?t:ooovri)- v] 3 ^ { 1 . 1377 3 + 0 . 3 1 2 69'R'' R1 1 -A 
l'-0.16254-" (41 + 1.''L: !l,/??.)‘-TE'’ 1 '-TFa 2+VV /-S/4 5 . 

AFA= ( 5,6 2C3S-''-( Ri + 1 . ' D/3 2. i 
141E ‘2-"Tl-:'Trl 2-'-\ ] l-rTE’'l-'TR12"'V12 
; Fc= ( (1. 127T5„-ra-rO.Ci 2CR9 ) '-C/:2. ) 

•’+TEf'2-Tl <7^12 ''-721 t-TE 'l'aR12«-V22 
ArC=(r7l-^ ( 1.13-’7a'';-pa + 'J.S126R‘^+Rl ) vD/3 2. ) 

1 + TEU2 -Tl->;-TPl2''-V31+TE*‘l-''-"'l-r EY'^V32 
A3A= ( (1. 1377S''-Rl-.-0.3:26P9 ) ; D/3 2, ) 

1 +TE* 2':Tl4a''’12x VVll-r-TE'‘l+T’'’12-';VV12 
-vGD= ( ( 1,787 936 '^Ri - .270000) "0/32 . ) 

1 +TE',2-:-T] \TR12-«-VV21-rTE '1MP12 + VV2Z 
AGC= (ra- t0.162 5 4C 'lRl+1. ) )'-D/32.) 

1 +TE,'12-!-T14 9 1-;: Vv'31-rTE U 'fR 1 x--'' E' 1 -« VV32 

AhA= ( ( i . 13 7 /e + 0.3 1269^'K-RI ) :-D/32 . ) 

1 +TEM2-/-Tl-'-TF12xvVVll + TE‘'a '.4T9 12;:VVV12 
mHB= (0.1625‘+0 MRl + l. )-':D/32. ) 

: +TEf 2 -;v-T‘'-;:TR12-::vVV21+TEN'1-'4TR124vVV 22 
AHC=(9 1^ (1.787 936 + 0.2 70'''00;'-pl)"'D/32. ) 
i + T E ‘ , 2 -; T 1 4 R 1 4: E M - V VV 3 1 - T E V 1 R 1 E ' "'C- V VV 3 2 
C 

C EVALUATIOF; OF VALJE3 OF A,3,C 

r 

PP =-aF' 4 ( MoL”-rtHC— AGC4 ah 3 ) + ^^0-!:- ( AFi^vAhC— AH34-AFC ) 

--A.t'' i AF2-4ASC-AGB“-ARC ) 

A H A 1- A H -'- « G L ) -AG'-":' ( ~AF AHC+AH :: \FC) 
l-'-AHA: (-AF-:>ASC + AC''-''FC) 

-AH4cAGo-f-AG*Arib ) -AGA^c ( -AHi AFB+AF* AH3 ) 

1 I AhA-;: (-AG-X-AFB+AF+ACB ) 

A J = /aF A-K- { ACB+AHC ~AHb+ACC ) + AGA'' ( AriB <!AFC~mF3-::-AHC ) 

1 + AMA ; ( AFB' AGC-AGE-'4AFC ) 



c 

c 

c 
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11 

20 

21 
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1 
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3 ' 

4 

5 

6 


3 

Q 


fp=p°/aj 

2Q=C3/AJ 
PP=RR/AJ 
A=A+FP 
3 — B + G Q 
C=C+RR 

PRlM18,,'.,i1,A ,B ,C 

F0P.'mT(31> ,I^s11X,I?,1(jX,E13.°,1CX,l15.C,1CX>E15.8} 
IF(r-ujSCP"').LT»l. E-6 Ab's ( 0'^ ) .LT • 1 • E-£ .Ar:D.AbS(Pci) .LT.l.E-6) 

iGOTOZo 
CONTIRJE 
PRliUZl 

For -"AT ( /?5>' »5HTHETA » 1 PX , 6MT0P0'JE , 2 3X , SHTJB ( XAX ) / ) 

CT = ‘-*-»''Rl'''(5» '■A + 3-t-Rl■5^C)/A5, 

DT0J = 2, "ra-' { A^5 ) 

RRlf T22,,1,DT,DT0J 

F 3k I',. A I (35X»I1G»15X»F1G«4»15X)F1G»4) 

CO'.TIk JE 
C0ImTIi'4UE 
STOP 
END 


CUbROUTINi- COiiFjTtiS VALuC CF 


J, V I L. b K H I w r 1 


GRID FONTS 


oOu nOl, I I HE oJwl ( A > . » > 3 , C , ,,2 5 >.3 5 aA 5 X5 j Y2 , Y3 , YA 5 Y5 > EM , R1 ,F7 » F8 s 
lF 9 ,Fio,rll,Fl 2 !.Fl 3 ,Fl 4 ,Fl 5 ,Fl 6 >Fl 7 >Fl 8 ,Fl 9 ,F 20 !.F 21 ,F 22 ,F 23 ,F 2 A> 


2FZS>F26»r<^~’) 

Fl = Rl:tX2\-VA-v( A+5-' X5 + R1 x C ^ 7 2 ) 2 +Y 2 -^ X aa ( /J + X2 + R 1 C- '^5 ) x 2 

IF(Fl)l,lC»Z 

PRINT2,^1 

FORMAT (lx, ElO. A) 

STOP 

IF(moS(X-1. ) .LT .1.E-06.AnC.ASS{ Y~1 • > . LT . 1 . E-C6 ) OCTOl 0 

IF { EM-1. ) A,5, A 

F2 = F1«^! ( EiM-1. ) 

oOTOo 

F2 = X , 


IF(EN-2.)7,8,7 
r3 = Fl"'f-“ { EY-2. ) 


C0T09 
F3 = l. 

FA=R1-.5'X2 - YA^H A + E-5:X5 + P1'.C*Y2 ) 
F5 = Y2:rX^-' ( A+B*X2 + R 1^-' 5 ) 
F6^FA+F5 
FT = F2 x-f6 

F6=F2*(FA^X5+F5*X2) 

F5 = F2-:'- ( FA-kyZ+FSts- f5 ) 

-10 = F 3 xFS"--i 2 

F11 = F2-c(K1-xx2--YA + y2^XA) 

F 12 = F 3 -xF 64 ( fA--X 3 + F 5 *X 2 ) 

F 13 = F 2 --^(R 1 -;>X 2 ^' Tk-xXS + YZ-tXAvXZ ) 
FlA = F3-;F6-:MFA-«Y2 + F5-!rY5 ) 

F15 = F2» {R1-;^X2-x Y4-W2+Y2«XAx rs ) 
F1S = F3X'Fo'> (FA^X5-rF5^!X2) 



1'' 


F 17 = F2-;i- ( R 1-.; X2 X '( X5 + '^2 ;'-y>^rX2 ) 

1 S = r 2 "- ( F V 5 4 F 5 X y 2 ) - 2 
F IS = - 2* ( R1 X X2*Y '4 xyS ''-X5-t-Y2 vy^^x2 x- - 2 ) 
F20='"3x(FSxY24-'"5xv5)';(p4<x54.p5-''2' 
-21 = F2i- ( Y4 X5XY24-V2 ''4--cY5 x-X2 ) 


F22 = F3-';-{ i'24F5 " /5 ) ~'’FS 

P23 = F2 - 1 Rl-;A2-'Y'+';Y2-r ,'2':xA ; V5 ) 

F2F=F2 ; ( FA 'X54F5 -'Xl ) - ( - “ V2 + F 5 / ''5 ) 

r2;-F2-x-(Rl''\2^' ''r' Y2^ X5-rY2AX4 

F26 = F2:-(FA:' 2-rF5. /S ) x''2 
F2-'^F2--(,s:--X2-->'^^ XY2-> Y24-YZXXA-X ,5-x , 5 ) 
R2TURR 


-7 = 0.0 
FS = o.O 
F9 = 0.0 
- 10 = 0.0 
F11 = 0.0 
F12=0.0 
F13=0,0 
F14=0.0 
F15=0.0 
FI. 5 = 0.0 
FI 7-0 ,0 
FlS=0.n 
F19=0.0 
F 2 0 = 0 • 0 
F21=0.0 
FZ2 = ri,p 

r2 3 = C • 0 
F24=0.0 
F 2 S = 0 . 0 
F26 = 0 .n 
F27=0.0 
return 


El'JD 


RTPY 



